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ABSTRACT 

The  purpose  of  the  Phase  11  program  entitled  "Single  Crystal  Growth  Optimization  of 
Magnesium-Doped  Lithium  Niobate"  was  to  optimize  the  growth  of  Mg0:LiNb03.  The 
first  step  in  this  direction  consisted  of  developing  means  and  methods  to  accurately 
measure  absorption  losses  (scattering  loss  methods  were  previously  developed  under  a 
Phase  1  contract).  It  was  known  at  the  outset  that  very  sensitive  instrumentation  would 
be  required  to  allow  measurements  on  small  samples.  The  results  of  the  measurements 
of  both  lypes  of  photon  loss  would  then  be  used  to  optimize  the  molar  concentration  of 
MgO  in  UNbOj  to  achieve  crystals  with  low  loss  as  well  as  high  photorefractive  damage 
threshold. 

Highly  sensitive  scattering  and  absorption  facilities  were  developed  and  used  in  the 
process  of  optimizing  crystal  growth.  The  ability  to  accurately  and  routinely  measure 
these  parameters  has  allowed  us  to  understand  variability  in  the  growth  process, 
determine  growth  parameters,  and  -most  importantly-  to  observe  that  specific 
propagation  and  polarization  directions  are  significantly  less  lossy  than  others. 

Studies  were  performed  on  the  growth,  poling,  and  resulting  properties  of  MgO-doped 
lithium  niobate.  Both  photorefractive  damage  threshold  and  NLO  phase-matching 
temperature  were  observed  to  vary  with  MgO  concentration.  For  one  crystal  doped  with 
a  high  MgO  concentration,  it  was  impossible  to  induce  measurable  photorefractive 
damage. 
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PROGRAM  OVERVIEW 


1.1  ABSTRACT 

The  purpose  of  the  Phase  II  program  entitled  "Single  Crystal  Growth  Optimization  of 
Magnesium-Doped  Lithium  Niobatc”  was  to  optimize  the  growth  of  MgOiLiNbOj.  The 
first  step  in  this  direction  consisted  of  developing  means  and  methods  to  accurately 
measure  absorption  losses  (scat'erirg  loss  methods  were  previously  developed  under  a 
Phase  I  contract).  It  was  known  a:  the  o’ltset  that  very  sensitive  instrumentation  would 
be  required  to  allow  measurements  on  small  samples.  The  results  of  the  measurements 
of  both  types  of  photon  loss  would  then  be  used  to  optimize  the  molar  concentration  of 
MgO  in  lJNb03  to  achieve  crystals  with  low  loss  as  well  as  high  photorefractive  damage 
threshold. 

Highly  sensitive  scattering  and  absorption  facilities  were  developed  and  used  in  the 
process  of  optimizing  crystal  growth.  The  ability  to  accurately  and  routinely  measure 
these  parameters  has  allowed  us  to  understand  variab*’  i;  -ii  the  growth  process, 
determine  growth  parameters,  and  -most  importantly-  to  observe  that  specific 
propagation  and  polarization  direction  5  are  significantly  less  lossy  than  others. 

Studies  were  performed  on  the  growth,  poling,  and  resulting  properties  of  MgO-doped 
lithium  niobate.  Both  photorefractive  damage  threshold  and  NLO  phase-matching 
temperature  were  observed  to  vai>'  with  MgO  concentration.  For  one  crystal  doped  with 
a  high  MgO  concentration,  it  was  impossible  to  induce  measurable  photorefractive 
damage. 
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BACKGROUND 


Single  crystals  of  high  optical  quality  have  an  increasingly  important  role  in  electro-optics 
as  these  devices  become  widespread  in  military  and  commercial  systems.  As  a 
consequence,  more  exacting  specifications  are  being  placed  on  these  components  in 
order  to  achieve  optimum  results. 

lithium  niobate  (liNbOj)  is  presently  used  for  electro-optic  Q-switches,  harmonic 
generators,  acousto-optic  devices,  transducers,  mcxlulators,  optical  parametric  oscillators 
and  amplifiers  and  optically  bistable  devices.  It  is  one  of  the  key  elements  in  Nd:YAG 
lasers  which  are  in  rangefinders  and  target  designators  used  by  the  Army.  The  optical 
quality  of  the  material  is  of  critical  importance  as  it  determines  the  reliability  and 
performance  of  these  devices;  at  present  most  failures  occur  within  this  optical  element. 
One  critical  parameter  is  the  optical  loss  caused  by  the  material  itself,  which  can  be  due 
to  either  scattering  or  absorption.  To  optimize  crystal  growth  it  becomes  essential  to 
precisely  measure  these  two  losses  in  order  to  correlate  them  to  growth  conditions. 
However,  before  the  present  program  was  in  place,  there  was  no  facility  available  to 
reliably  and  routinely  measure  the  low  levels  of  optical  loss  in  these  crystals. 

Through  the  Czochralski  process,  crystals  as  large  as  150  mm  in  diameter  can  be 
reproducibly  grown.  High-volume  crystal  fabrication  techniques  have  been  developed 
enabling  high-yield  production  to  demanding  tolerances.  The  use  of  lithium  niobate  in 
optical  applications  in  the  near-IR  and  visible  regions  of  the  spectrum  has  historically 
been  limited  due  to  a  phenomenon  known  as  the  photorefractive  effect  (PRE).  In  PRE, 
the  laser  beam  itself  causes  the  uncontrolled  redistribution  of  space  charges  in  a  crystal. 
Through  the  crystal’s  electro-optic  effect,  this  redistribution  gives  rise  to  an  uncontrolled 
refractive  index  variation  that  in  turn  distorts  the  laser  beam.  Although  such  PRE 
damage  generally  does  not  lead  to  catastrophic  crystal  failure  and  can  in  fact  be  removed 
through  armeaiing,  its  uncontrollabUity  in  most  cases  renders  a  subject  material 
unsuitable  for  practical  device  application. 

The  doping  of  UNbOj  with  several  mol%  of  magnesium  oxide  (MgO)  was  first  reported 
in  1980  *  as  a  means  to  increase  the  material’s  PRE  damage  threshold.  Early  work 
focused  on  the  effect  of  MgO  doping  concentrations  and  identified  crystal  homogeneity 
problems  associated  with  a  non-unity  Mg  distribution  coefficient. 
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13  OBJECTIVES 


Iliere  were  three  major  objectives  in  this  Phase  II  program: 

1)  Develop  and  demonstrate  the  teclmoio^  to  routinely  measure  the  absorption 
losses  in  transparent  materials  with  absorption  coefficients  as  low  as  IC*  m■^ 

2)  j^ply  the  technology  developed  in  the  first  objective  to  measure  the  optical  loss 
in  Mg0:LiNb03  single  crystals. 

3)  To  optimize  the  growth  of  MgOrLiNbOj  for  minimum  optical  loss  and  maximized 
phoiorefractive  damage  threshold. 


1.4  ACCOMPLISHMENTS 

These  three  tasks  have  been  successfiilly  completed,  as  is  shown  in  the  body  of  this 

report.  As  a  result,  we  have: 

o  Developed  hardware  and  methods  to  measure  loss  in  crystals  reliably  and 

routinely.  The  accuracy  can  be  as  low  as  5  x  10'^  for  samples  that  absorb 
about  1%  at  modest  power  levels.  Extension  to  accuracies  of  10'^  is  possible  with 
higher  laser  power  and  larger  samples,  due  to  the  high  resolution  of  the  apparatus 
(4  X  10-^  °K). 

The  equipment  and  methods  are  general  enough  that  they  can  be  applied  to  any 
transparent  material  and  so  could  be  used  to  measure  the  small  scattering  and 
absorption  losses  in  a  wide  variety  of  materials.  This  information  could  then  be 
used  to  further  improve  these  materials  as  well,  and  should  lead  to  an  important 
reduction  in  the  development  time. 

o  Optimized  the  growth  parameters  for  MgOtLiNbOj  using  the  Czochralski  method, 
thus  making  it  possible  to  obtain  single  crystals  that  have  optimal  absorption  and 
PRE  characteristics. 

o  Identified  directions  of  propagation  and  polarization  that  cause  the  least  amount 
of  absorption,  and  at  the  same  time  determined  that  scattering  is  a  minor 
component  of  optical  loss  in  MgOtLiNbOj. 
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2  DEVELOPMENT  OF  OPTICAL  LOSS  MEASUREMENT  CAPABILITY 


Optical  leases,  as  mentioned  earlier,  can  be  divided  into  scattering  and  absorption.  A 
method  was  developed  during  the  Phase  I  of  this  program  to  measure  scattering  loss. 
During  the  Phase  U  this  method  was  further  developed  for  routine  testing,  although  no 
substantial  improvements  were  required. 


2.1  SCATTERING  FACILITY 

Our  scattering  loss  measurement  system  is  similar  to  a  system  recently  completed  at 
Stanford  University.^  The  system’s  major  components  are  an  integrating  sphere,  a 
chopper,  and  a  lock-in  amplifier,  as  shown  schematically  in  figure  1.  The  laser  in  the 


diagram  be  any  of  several  types.  We  have  used  He-Ne,  argon  ion,  and  Nd:YAG 
i£^ers  during  the  course  of  our  investigation.  The  laser  beam  is  sent  through  a  spatial 
filter  in  order  to  insure  a  clean  Gaussian  beam.  This  is  necessary  so  that  the  beam  will 
be  well  defined,  can  be  focused  to  a  well-defined  spot  and  will  be  as  noise-free  as 
possible.  The  laser  beam  is  then  sent  through  the  optical  chopper  which  also  sends  a 
reference  signal  to  the  lock-in  detector.  A  small  fraction  of  the  beam  is  picked  off  and 
sent  to  a  laser  beam  monitor  detector.  The  beam  is  then  focused  into  the  sample  aystai, 
while  apertures  are  placed  in  the  beam  to  block  out  any  light  that  may  have  been 
scattered  by  any  of  the  optical  elements  in  the  beam.  TTie  laser  beam  then  passes 
through  the  sample  crystal.  The  portion  of  the  beam  that  is  scattered  by  the  crystal  is 
captured  and  reflected  by  the  integrating  sphere  and  is  eventually  detected  by  a  baffled 
semiconductor  detector.  We  have  used  both  silicon  PIN  diodes  and  InGaAsP  detectors. 
The  output  from  the  detector  is  then  amplified  and  sent  to  the  lock-in  amplifier.  With 
this  s^tem,  depending  on  the  laser  used  and  the  laser  output  power,  we  have  been  able 
to  detect  a  minimum  signal  of  from  125  nanovolts  to  1.26  microvolts  with  laser  input 
powers  of  from  170  /tW  to  1.6  mW.  A  photograph  of  the  apparatus  is  shown  in  figure  2 
and  a  photograph  of  the  opened  integrating  sphere  is  shown  in  figure  3. 


F^ure  2  Photograph  of  the  scattering  facility. 
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2.1.1  Scattering  measurements  ••  methods. 


The  measurement  technique  involves  placing  the  material  sample  within  the  integrating 
sphere  and  focmsing  the  laser  radiation  through  the  sample.  If  both  sample  end  faces 
are  not  allowed  to  reflect  light  directly  into  the  integrating  sphere,  only  light  scattered 
within  the  sample  bulk  is  measured  by  the  detector. 


Figure  5.  Photoguph  of  the  integrating  sphere  for  scattering 
measurements. 


Some  early  experimentation  with  index  matching  oils  on  the  faces  to  eliminate  back 
reflections  into  the  integrating  sphere  indicated  that  the  need  for  "zero’’  scattering  light 
guides  does  not  warrant  using  these  oils.  The  difficulties  encountered  in  handling  the 
samples  and  keeping  the  integrating  sphere  free  from  oil  do  not  compensate  the  possible 
gains  in  measurement  accuracy.  We  chose  to  try  to  separate  that  part  of  the  detected 
signal  due  to  this  backscattered  light  from  the  total  signal  by  measuring  the  scattered 
signal  for  a  central  region  of  the  crystal  within  the  integrating  sphere.  Internal  baffles 
and  highly  reflective  internal  surfaces  direct  the  scattered  light  toward  a  detector  which  is 
mounted  on  the  integrating  sphere.  The  signal  measured  on  this  detector  gives  the 
relative  magnitude  of  the  scattered  light. 

Blocking  the  exit  port  and  measuring  the  detector’s  signal  gives  a  relative  measure  of  the 
total  power  into  the  sample.  Comparison  of  the  signal  obtained  with  the  sample  in  and 
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out  of  the  sphere  gives  an  absolute  determination  of  the  optical  scattering  losses. 
Chopping  the  laser  allows  phase-sensitive  detection  which  enhances  the  signal-to-noise 
ratio  at  low  scattering  levels.  With  the  appropriate  detector  and  reflective  coatings 
within  the  integrating  sphere  any  CW  laser  can  be  used  in  this  test  apparatus. 

We  should  point  out  that  this  apparatus  is  very  flexible  and  can  be  used  to  measure  the 
scattering  losses  of  any  transparent  material. 


2.12  Scattering  measurements  ~  data 

Because  the  samples  showed  a  very  small  scattering  coefficient  (about  0.005  1/m),  we 
opted  to  use  a  He*Ne  laser  beam  rather  than  a  YAG  laser  due  to  the  higher  cross 
section  at  shorter  wavelengths.  The  extreme  sensitivity  of  the  sensor,  on  the  other  hand, 
made  it  more  convenient  to  use  this  low  power  source  instead  of  a  10  W  Argon  ion 
laser.  All  the  measurements  reported  here  were  performed  with  a  5  mW  He-Ne  laser 
(632.8  nm)  beam  focused  to  approximately  50  in  diameter. 

Figure  4  shows  the  atomic  arrangement  of  UNbOj,  and  the  measured  scattering 
coefficients  are  summarized  in  table  1.  The  measurements  of  the  first  z-grown  crystals 
available,  Z-1  through  Z-4,  were  found  later  to  be  in  error  due  to  a  saturation  of  the 
amplifier.  At  that  point  a  new  formulation  was  announced,  whereby  a  more  closely 
congruent  mixture  of  Li  to  Nb  would  be  used,  and  for  this  reason  the  measurements 
were  not  redone  in  these  samples.  Having  corrected  the  saturation  problem,  we 
performed  the  measurements  on  crystals  prepared  from  boules  Z-5  on. 

No  clear  differences  could  be  noticed  between  Z-,  Y-  or  X-grown  crystals,  except  for  the 
obvious  defects  that  could  be  observed  in  the  boules  (Y-4  and  Y-5,  for  instance).  There 
is  a  clear  difference,  however,  in  the  scattering  observed  for  Z-propagation  compared  to 
any  of  the  other  two  directions.  Although  a  slight  difference  is  observed  in  the 
absorption  coefficients  (see  below)  with  different  polarizations,  this  was  not  observed  for 
scattering.  Contrary'  to  expectations,  striae  were  not  found  to  have  a  major  effect  on 
scattering,  suggesting  that  the  index  modulation  is  very  small  compared  to  other  causes. 
Figure  5  contains  the  data  from  table  1,  having  grouped  all  polarization  and  growth  axis 
together,  and  differentiating  only  on  the  basis  of  propagation  directions:  either  Z-  or 
noii-Z  axis  propagation. 


10 


Schematic  atomic  arrangement  in  the  basal  (0001)  plane  of  lithium  niobate. 
The  oxygen  subiat^^c^  b  represented  by  lines  connecting  nearest-neighbor 
oxygen  centers;  the  two  layers  of  oxygen  atoms  sandwich  the  cations  into 
octahedral  sites.  N^ne  that  there  exist  three  symmetrically  equivalent 
sets  of  Cartesian  x-  and  y-axes;  only  one  is  shown. 


Fiffjure  4.  Cartesian  axis  definition  for  LiNbO^ 


Scottcring  cMfRcwnt  [t/ni| 


Table  1.  Summary  of  scattering  coefficients  [l/m]  at  632.8  nm. 


PROPAGATION  > 

■z 

Ss*Si=?a 

XOR  Y 

POLARIZATION  > 

XORV 

YORX 

Z 

•ANY 

BOULE  [Lip] 

[MgO] 

AVG 

AVG 

AVG 

AVG 

Z-5 

iS 

4.65 

0.027 

0.006 

0.008 

0.007 

Z-6 

48.38 

4.70 

0.032 

0.017 

0.021 

0.019 

Z-7 

48.38 

4.80 

0.017 

0.005 

0.006 

0.005 

Z-8 

48.38 

4.90 

0.025 

0.006 

0.005 

0.005 

X-3 

48.38 

3.00 

0.019 

0.006 

0.008 

0.007 

Y-1 

48.38 

4.70 

0.019 

0.006 

0.009 

0.008 

Y.2 

48.38 

4.80 

0.011 

0.009 

0.009 

0.009 

Y-3 

48.38 

3.00 

0.017 

0.008 

0.008 

0.008 

Y-4 

48.38 

4.90 

0.105 

0.061 

0.049 

0.055 

Y-5 

48.38 

5.10 

0.149 

0.139 

0.124 

0.132 

Avg  indicates  average  over  all  samples  in  each  boule. 


22  CALORIMETRIC  PHOTO-ABSORPTION  FACILITY 


22.1  Calorimetiy  of  photo-absorption 

The  flow  of  energy  into  and  out  of  a  crystal  cell  is  schematically  shown  in  figure  6.  The 


Figure  6  Energy  flow  diagram  for  photoabsorption  in  a  crystal. 


index  of  refraction  n  determines  the  amount  of  energy  that  enters  the  crystal,  Pj^,  and  the 
fraction  converted  into  heat  (inelastic  scattering)  is  k  Another  fraction  is  elastically 
scattered,  leaving  the  cell.  As  a  consequence  of  these  two  losses,  the  beam  intensity 
decreases  exponentially  along  its  path: 


I(x)  =  lo  e 


-  (A  +  a  )  X 


where  the  constants  A.  and  a  correspond  to  absorption  and  scattering,  respectively.  The 
energy  deposited  in  the  crystal  by  absorption  raises  its  temperature  as  a  function  of  time 
The  absorption  coefficient  k,  is  given  by  the  path  length  L  and  the  converted  fraction  k 
as 

k  =  (l/L)ln[(T-RK  )/(T-k)] 
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where  T  and  R  are  the  transmission  and  reflection  coefficients  for  the  crystal  face.  The 
appendices  contain  details  of  this  derivation  and  the  means  to  convert  the  time 
dependence  of  temperature  into  x. 


Measuring  the  absorption  coefficient  was  accomplished  with  a  very  sensitive  calorimetric 
photo  absorption  apparatus,  described  schematiiklly  in  figure  7. 


Figure  7.  Schematics  of  the  calorimetric  absorption  spectroscopy  setup. 


A  vacuum  chamber,  in  conjunction  with  a  turbomolecular  vacuum  pump,  create  an 
environment  free  of  convection  and  air  conduction  losses  for  the  samples,  by  sustaining 
them  at  pressures  below  Torr.  Light  generated  by  a  laser  is  stabilized  by  a  Laser 
Power  Controller  (LPC),  to  insure  that  the  delivered  power  is  constant  throughout  the 
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active  absorption  period.  The  beam  is  then  collimated  and  focused  into  the  samples 
(numbered  1  to  4)  and  its  intensity  is  monitored  by  the  laser  power  meter. 

The  vacuum  chamber  is  made  of  stairJess  steel  with  entrance  and  exit  windows,  a  top- 
view  window  to  facilitate  sample  location,  a  mechanical  feedthrough  to  move  the  samples 
in  and  out  of  the  beam,  and  electrical  feedthroughs  for  temperature  monitoring. 

The  sample  holder  (see  figure  8)  consists  of  a  frame  on  which  four  samples  can  be 
supported  by  a  grid  of  60  urn.  nylon  fibers.  The  limited  contact  surface  not  only  reduces 
the  loss,  but  also  reduces  variations  among  samples  of  different  mass. 


1800 


1 


Figure  8.  Frame  used  to  hold  samples  in  absorption  measurements.  Dimensions  in  inches. 

To  prevent  the  beam  from  touching  the  other  walls  of  the  cell  or  the  thermistor  itself,  it 
is  important  to  ensure  that  the  beam  diameter  is  safely  contained  within  the  sample, 
including  at  least  two  of  the  Fresnel  reflections.  This  was  achieved  by  having  a  beam 
diameter  of  less  than  500  pm,  while  the  crystal  dimensions  were  at  least  6x7  mm.  The 
inclusion  of  the  Fresnel  reflections  was  insured  during  sample  alignment  using  a  He-Ne 
pilot  beam. 

The  temperature  of  each  crystal  sample  is  monitored  by  low-mass  thermistors,  typically 
60  mg,  attached  to  the  middle  of  one  of  the  sides  with  Tra-Con  2902  Ag  epoxy  (for  good 
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Figure  9.  Photograph  of  calorimetric  absorption  facility. 


thermal  conductivity).  As  the  mass  of  the  crystals  used  was  between  4  and  8  grams,  a 
small  correction  (always  below  0.6%)  was  made  for  the  thermistor /epoxy  mass 
contribution  (see  Appendix  A  for  details).  Tlie  thermistors  chosen,  YSI  44037  and 
Alpha  Thermistor  Inc.,  part  #  13A6001-D3,  have  a  nominal  resistance  of  6  Kn  at  25  °C 
(to  minimize  lead  resistance  effects),  and  high  SR/ST  sensitivity.  Long,  thin  leads  (.004" 
diam.)  attach  the  thermistors  to  the  apparatus  to  minimize  loss  through  this  thermal 
cormection. 

The  resistance  is  mea.sured  (accurate  within  the  equivalent  of  4x  10'^  °C)  with  a  scanning 
ohmmeter  (Azonix  ScannerPlus  equipped  with  thermistor  cards)  that  measures  voltage 
generated  by  a  lO^uA  DC  current.  It  then  converts  the  measured  value  to  ASCII  code 
for  transmission  via  RS-232  to  a  PC.  The  Azonix  scanning  ohmmeter  is  capable  of 
sampling  up  to  four  channels  every  two  seconds. 

Two  light  sources  were  employed.  Initially,  an  Argon-ion  laser  tuned  to  476,  488  and 
514.5  nm  was  used  and  later  a  YAG  laser  was  employed  for  the  bulk  of  measurement  as 
absorption  at  1064  nm  is  more  relevant  to  LiNbOj  applications.  The  beam  was  normally 
polarized  (95%)  along  a  horizontal  plane,  and  sensitivity  to  polarization  was  studied  by 
rotating  the  polarization  with  a  A.  /2  plate. 

Laser  power  (in  the  range  of  200  -  500  mW)  was  monitored  with  an  Ophir  Optics  DG 


Fiffire  10.  Photograph  of  the  absorption  facility  vacuum  chamber. 


laser  power  meter.  Its  stated  accuracy  at  the  operating  range  used  is  5  mW,  equipped 
with  A-D-Serial  c’-  'luitry  for  logging  the  data  on  the  PC. 


2.2.2  Measurements 

Six  parameters  are  monitored  by  the  computer:  time  (seconds),  laser  power  (mW),  and 
resistances  of  all  four  thermistors  (Ohms).  The  data  acquisition  program  performs 
measurements  about  every  500  ms,  and  records  a  measurement  (writes  it  into  a  data  file) 
if  any  of  the  above  parameters  differs  from  the  prior  recorded  value  by  more  than  a 
presettable  threshold.  This  allows  faster  sampling  when  the  temperatures  change  rapidly, 
or  when  power  variations  are  detected  (turning  laser  on  or  off),  or  at  maximum  time 
intervals,  while  keeping  the  amount  of  data  to  a  manageable  size. 

The  program  computes  temperatures  from  the  resistance  values  using  the  Steinhart  and 
Hai  t  equatioi.’  and  coefficients  obtained  from  the  calibration  curves  provided  by  the 
manufacturers  (YSI  and  Alpha  Thermistor  Inc.).  These  four  temperature  values  are 
appended  to  the  six  raw  data  points  for  each  observation. 
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The  data  analysis  uses  the  statistical  methods  outlined  in  appendix  A,  based  on  the 
following  experimental  conditions: 

•  Laser  power  is  kept  constant  during  the  heating  period  to  within  1% 

•  Vacuum  is  kept  below  lO"^  Torr  throughout  heating  and  cooling 

•  Heating  is  applied  for  a  time  at  least  as  long  as  the  cooling  time  constant. 

Preliminary  runs  were  made  using  two  thermistors  attached  to  a  single  crystal,  one  in  the 
center  and  one  near  the  entrance  point,  as  shown  in  figure  11. 


Fiffire  11.  Two-thennistor  configuration  used  during  system  checkout. 
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The  resulting  thermal  measurement  illustrates  the  level  of  sensitivity  of  the 
instrumentation.  We  measured  four  temperature  values  (two  per  crystal)  on  two  crystals 
of  the  same  dimensions.  The  first  crystal,  monitored  by  thermistors  1  and  2,  was  not 
exposed  to  the  beam.  The  second  crystal  (monitored  by  thermistors  3  and  4)  was 
exposed  to  a  13  W  (CW)  beam  at  1064  nm.  As  shown  in  figure  12,  the  beam  was  on 
from  t  *  20  to  t  =  1400  seconds,  and  was  quite  stable. 


Figure  J2  Typical  measurement  of  laser  power. 


Figure  13.  Raw  temperatures  for  all  four  thennistors  depicted  in  figure  11. 
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The  raw  temperatures  shown  in  figure  13  indicate  a  slight  drift  throughout  the  chamber, 
as  indicated  by  the  reference  thermistors  1  and  2. 


F^ure  14.  Temperatures  near  the  front  surface  (T3-T1)  and  at  the  center  of  the  crystal  (T4-T2). 


Figure  15.  Temperature  gradient:  front  face  to  middle  of  crystal  sample. 


When  the  reference  temperatures  are  subtracted,  the  temperature  changes  due  to  the 
beam  can  be  observed,  as  shown  in  figure  14.  Both  curves  (T3-T1)  and  (T4-T2)  match 
quite  closely  except  for  a  small  offset  due  to  variations  intrinsic  to  the  sensors.  Tins 
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offset  can  be  observed  under  nc  power,  averaged  over  some  time,  and  subtracted  from 
each  thermistor’s  readings. 

The  difference  (T3-T4)-(T1-T2)  is  plotted  in  an  enlarged  scale  in  figure  15,  where  a 
slight  temperature  gradient  along  the  crystal  can  be  observed.  The  curve  shows  that  the 
end  of  the  rod  (T3)  warms  up  more  rapidly  than  the  bulk  (T4)  at  the  beginning,  probably 
due  to  surface  imperfections,  until  equilibrium  is  reached  and  thereafter  the  end  stays 
cooler  than  the  bulk  (more  access  to  radiative  surface).  The  effect  is  reversed  when  the 
laser  is  turned  off. 

Fi^re  16  shows  a  typical  computer  printout  for  a  single  measurement  of  absorption 
coefficient  The  values  used  for  all  the  calculations  are  present  in  each  printout  as  well 
as  graphical  information  on  the  exponential  fits. 


223  Uncertainties 

Several  uncertainties  have  a  role  in  the  measurement  of  the  absorption  coefficient  X,  and 
not  all  contribute  in  a  predictable  way  to  the  final  accuracy  of  the  measurements.  Direct 
measurement  errors  come  from  the  following  sources: 

Pq  •  Incident  power.  Known  to  1%  relative,  3%  absolute. 

Cg  -  Crystal  heat  capacity.  Not  known  for  MgO  doped  LiNbOj,  and  assumed  the  same 
for  all  crystals  to  be  0.654  tJ/g-“C],  the  value  for  undoped  liNb03. 

L  ’  Light  path  length.  Typically  25  mm,  known  ±  5  /im. 
n  -  Ciy'stal  index  of  refraction.  Assumed  for  all  c.ystals  to  be  2.25. 
m  -  Crystal  mass.  Typically  4  -  6  g,  known  within  ±  2mg. 

Other  aspects  contribute  to  the  total  uncertainty,  such  as  cleanliness  of  the  crystal  and 
window  surfaces,  beam  purity  and  alignmeni,  and  drift  in  the  measuring  equipment.  As 
expected,  these  less  quantifiable  components  of  the  total  uncertainty  contribute  the  most. 
As  a  test  of  reproducibility,  the  same  sample  was  reinstalled  and  measured  throughout 
the  14  months  of  measurements,  yielding  a  value  that  oscillated  between  0.194  and  0.199 
[1/m].  We  believe  that  this  represents  the  experimental  (random)  uncertainty  in  the 
results,  estimated  for  all  samples  at  aoout  0.005  [1/m].  Systematic  errors  (for  example,  if 
the  assumed  values  for  C^  and  n  were  fennei  in  error)  affect  all  measurements  *-  and  can 
be  corrected  for,  should  new  values  appe^iT  more  accurate.  For  the  purpose  of 
optimizing  crystal  growth,  however,  systematic  errors  do  not  alter  the  conclusions  as  the 
dependency  of  X  is  monotonic. 
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Fi^ire  16.  Typical  computer  printout  for  absorption  measurements. 
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23  RESULTS 


Table  2  sho^  the  samples  obtained  for  each  grown  boule,  their  corresponding 
(impositions,  and  the  measured  values  of  absorption  coefficient  along  the  various  axis, 
differentiating  among  1)  propagation  along  the  Z  axis,  any  polarization,  2)  propagation 
perpendicular  to  the  Z  axis,  polarization  parallel  to  Z  axis,  and  3)  propagation 

2  Summary  of  absorption  coefficients  [1/m]  ai  1064  run. 


PROPAGATION  > 

as 'SI  a  .ss.:as 

xo 

POLARIZATION  > 

XORY 

YORK 

BOULE  [UfO] 

[MgO] 

AVG 

u 

AVG 

u 

AVG  U 

Z-1 

4.80 

0231 

6.6 

0.075 

73 

0.037  2.6 

Z.2 

4.70 

0.109 

0.052 

Z.3 

4.60 

0.155 

0.110 

52 

0.050  83 

Z4 

4S.60 

430 

0.179 

0.089 

0.046  3.8 

Z.5 

4.65 

0.133 

7.4 

0.057 

0.044 

Z-6 

4838 

4.70 

0.175 

23 

0.112 

0.053 

Z.7 

4838 

4.80 

0,150 

0.116 

Z-8 

48.38 

4.90 

0.127 

0.101 

X-3 

4838 

3.00 

0.158 

0.132 

Y-1 

4838 

4.70 

0,164 

0.131 

Y.2 

4838 

4.80 

0.266 

3.6 

0.246 

0.114 

Y-3 

48.38 

3.00 

0351 

8.6 

0.222  1.6 

Y-4 

48.38 

4.90 

0.124 

0.152 

Y-S 

48.38 

5.10 

0.547 

0.408 

Avg  indicates  average  over  all  samples  in  each  boule. 

U  ®  Uniformity  rating:  average  value  divided  by  observed  range,  n>3. 


perpendicular  to  the  Z  axis,  polarization  perpendicular  to  Z  axis. 

The  same  information  is  shown  in  figures  17  through  21.  The  two  LijO  concentrations 
of  48.60%  and  48.38%  in  the  case  of  Z-grown  crystals  are  plotted  in  figures  17  to  18. 

The  more  closely  congruent  composition  of  48.38%  Li^O  was  used  for  the  Y-grown 
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AesoRPTioN  cocrriccNT  z -crown 

U}0  46.C0  fneiX 


Figure  17.  Absorption  coefficient  for  Z-grown  boules  using  the  quasi-congruent  mixture  with 
48.60  mol%  Li20. 


ABSORPTION  COCPnCCNT  2-GROWN 
tt20  4t.5a  mors 


Figure  18  Absorption  coefficient  for  Z^grown  boules  using  the  newly  developed  mixture  with 
48.38  mol%  U2O. 

(  'stals,  shown  in  figures  19  and  20.  In  these  graphs  the  measurement  error  bars  are 
aj,  oximately  the  size  of  the  symbols. 
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ABSORPTION  COErrOCNT  2-GROWN 
U20  48.«0  A>0  4848  moR 


{  ■  I  PttO^AiCH  ^  QT»C»  ^  O  Ott^H/OTMP 


Figfve  19.  Absorption  coefficient  for  Z-ff’own  crystal  Data  for  48.60  and  48.38  mol% 
combined. 


absorption  cocrriccNT  Y-o:?owN 

LQO  48.38 


I  m  i  Mo^AQAtioH  iA  OTIC*  mn^/z  ot*^9/oiwt 


Figure  20.  Absorption  coefficient  for  Y-grown  boules  with  48.38  mol%  Li20.  Subset  plotted  in 
the  same  scale  as  figure  19  for  comparison. 
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Figure  21.  Absorption  coefficient  for  Y-grown  LiNbOj. 


Comparing  figures  17  and  18  does  not  indicate  that  the  different  Lij^  percentage  makes 
a  noticeable  difference.  Based  on  this  observation,  data  from  both  formulations  were 
combined  into  figure  19,  where  an  interesting  trend  can  be  seen.  It  appears  that  Z- 
propagation  causes  the  most  absorption,  whereas  Z-polarization  causes  the  least 
absorption.  This  is  surprising,  since  the  striae  are  perpendicular  to  the  Z>axis  and  it  was 
believed  that  they  were  a  major  cause  for  the  increased  absorption  in  MgOiLiNbOj.  It  is 
also  interesting  that  within  the  range  of  MgO  concentrations  between  4.5  and  5  mol% 
there  is  no  visible  increase  or  decrease  in  absorbance. 

A  subset  of  the  Y-grown  crystal  data  was  plotted  in  figure  20,  ai  the  same  scale  as  figure 
19,  for  comparison  purposes.  Although  the  data  are  scanty,  no  dramatic  difference 
appears  to  exist  between  Z-  and  Y-grown  crystals,  except  perhaps  in  the  fact  that  Y- 
grown  seemed  more  difficult  to  obtain  without  cracks  or  other  defects  (see  Table  1  and 
figure  5  for  scattering  results). 

Figure  21  contains  all  the  data  obtained  from  Y-grown  crystals,  in  which  some  extreme 
MgO  concentrations  (from  3%  to  5.1%)  were  attempted.  The  5.1%  data  (Y-5)  show 
increased  absorption,  although  it  is  not  clear  whether  this  is  due  to  the  increased 
scattering  that  this  boule  exhibits. 
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3  CRYSTAL  GROWTH  DEVELOPMENT 


The  doping  of  LiNbOj  with  several  mol%  of  magnesium  oxide  (MgO)  was  first  reported 
in  1980  ^  as  a  means  to  increase  the  material’s  PRE  damage  threshold.  Early  work 
focused  on  the  ehect  of  MgO  doping  concentrations  and  identified  crystal  homogeneity 
problems  associated  with  a  non-unity  Mg  distribution  coefficient.  All  of  the  MgO-doped 
lithium  niobate  crystals  reported  to  date  with  MgO  concentrations  in  the  vicinity  of  5 
moi%  have  been  found  to  contain  planar  striae,  perpendicular  to  the  growth  direction 
with  thickness  and  spacing  on  a  su^millimeter  scale.  These  striae  appear  to  be  due  to 
microscopic  fluctuations  in  the  MgO  concentration,  corresponding  to  periodic 
microscopic  instabilities  in  the  crystal  growth  and  dopant  incorporation  processes. 

Without  engaging  in  a  broad-based  research  effort  to  investigate  the  details  of  the 
time-dependent  dopant  incorporation  phenomena  occurring  at  the  gi’owth  interface,  it 
seems  unlikely  that  the  formation  of  such  striae  could  be  eliminated  or  even  significantly 
limited. 

For  certain  applications  of  interest,  optical  interactions  in  directions  parallel  to  the  striae 
can  result  in  undesirable  beam  scattering  effects.  Of  particular  interest  is  noncritical 
phasematching  to  achieve  frequency  doubling  of  a  l.(^  fim  Nd:YAG  laser  beam.  It  has 
been  shown  that  such  an  application  would  preferably  be  addressed  through  use  of  a 
doubling  crystal  with  its  striae  oriented  parallel  to  the  Z-direction.  Such  a  crystal  would, 
therefore,  require  its  growth  axis  to  lie  in  the  X-Y  plane. 

Due  to  the  non-unity  distribution  coefficient  of  the  Mg  dopant,  MgO-doped  UNbOj 
crystals  exhibit  both  short-range  and  long-range  composition  variations.  The  short-range 
variations  took  the  form  of  periodic  planar  striae  oriented  perpendicular  to  a  crystal’s 
gro'vth  direction  with  a  spacing  of  roughly  200  fim.  These  striae  result  in  severe  scatter 
losses  for  laser  propagation  in  directions  perpendicular  to  a  crystal’s  growth  axis.  On  the 
other  hand,  a  long-range  composition  variation  was  observed  to  result  in  a  slow  change 
in  a  crystal’s  birefringence  along  the  direction  of  its  growth  axis.  This  in  turn  was 
observed  to  limit  the  effective  usable  crystal  length  for  laser  propagation  in  the  direction 
parallel  to  the  crystal’s  growth  axis.  It  thus  appeared  that  the  incongruency  of  the  growth 
process  for  Mg0:LiNb03  might  entail  inherent  limitations  for  application-related 
material  quality. 

The  aim  of  the  present  work  was  to  study  the  growth,  poling,  and  resulting  properties  of 
Mg0:LiNb03.  Ferroelectric  poling  conditions  were  studied  and  refined.  Growth-induced 
striae  were  observed  for  all  crystal  orientations.  Achievable  freedom  from  grown-in 
strain  was  found  to  depend  strongly  on  growth  orientation.  Over  the  range  of  conditions 
studied,  optical  losses  were  found  to  not  depend  strongly  on  MgO  concentration  or 
growth  orientation.  Photorefractive  damage  threshold  was  observed  to  be  a  function  of 
temperature  and  laser  fluence.  Growth  parameters  investigated  included  MgO 
concentration  in  the  melt,  growth  orientation,  and  basic  Czochralski  growth  conditions. 
Ferroelectric  poling  conditions  were  studied  and  refined.  Crystals  were  analyzed  for 
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MgO  concentration,  strain  state,  IR  transmission,  scattering  and  absorption  loss,  PRE 
damage  threshold,  and  nonlinear  optical  (NLO)  phase-matching  properties.  The 
discontinuous  shift  in  the  IR  absorption  peak  associated  with  residual  OH'  appeared  to 
be  an  important  indicator  of  dramatic  optical  property  changes  in  the  material. 


3.1  CRYSTAL  GROWTH  RUNS 

Sixteen  growth  runs  were  performed,  all  by  the  Czochralski  technique.  The  resulting 
boules  averaged  roughly  80  mm  in  diameter  by  50  mm  in  length.  Basic  charge 
preparation  and  crystal  growth  procedures  were  essentially  the  same  as  those  used  in  the 
production  of  optical-grade  undoped  LiNb03.  MgO  powders  were  screened  through 
spark  source  mass  spectrometry  to  ensure  chemical  purity,  after  a  drying  procedure  was 
refined  for  this  material. 

Table  3  summarizes  the  crystal  growth  experiments.  The  third  letter  in  the  "Boule  #" 
designation  indicates  the  growth  axis.  Crystals  were  grown  along  the  x,  y,  and  z  axes. 

Two  different  Li/Nb  ratios  were  used  in  the  course  of  the  experiments.  The  48.60  mol% 
LijC^  composition  had  been  the  industry  standard  for  "congruent’*  undoped  LiNbOj  at  the 
time  of  this  project’s  initiation.  While  this  project  was  underway,  parallel  work  at  CTI 
identified  a  corrected  value  of  48.38  mol%  for  the  actual  congruent  composition  for 
undoped  LiNb03  This  modified  value  was  then  adopted  for  subsequent  Mg0:LiNb03 
work.  A  variety  of  MgO  concentrations  were  used  in  doping  the  crystal  growth  melts, 
from  3.0  to  5.1  mol%.  Several  ciystals  were  analyzed  for  MgO  concentration  using  DC 
plasma  spectroscopy. 

IR  transmission  scans  were  performed  on  each  boule  as  a  function  of  axial  position, 
using  a  Perkin-Elmer  Lambda-9  Spectrophotometer.  The  technique  involved  using 
masking  tape  to  define  3mm  square  windows  immediately  below  the  shoulder  of  a  boule 
and  immediately  above  the  foot. 

Ferroelectric  Curie  temperature  measurements  were  performed  on  samples  from  the  top 
and  bottom  of  each  boule.  The  measurement  apparatus,  based  on  dielectric  effects,  has 
been  described  previously  *  and  is  routinely  used  at  CTI  to  investigate  compositional  and 
poling  issues  in  ferroelectric  crystals. 

Czochralski  growth  conditions  proprietary  to  CTI  were  initially  based  on  those  used  in 
production  of  optical-grade  undoped  LiNbOj  and  were  then  adjusted  in  an  iterative 
manner.  Immediately  following  growth,  each  boule  was  inspected  for  macroscopic  crystal 
quality  using  .standard  production  procedures  for  undoped  LiNbOj.  In  Table  3,  "Stiain' 
is  rated  on  a  qualitative  scale  of  A,  B,  C,  D.  F  with  "A"  representing  typical  undoped 
LiNb03  suitable  for  0-switch  application  per  CTI  specifications,  "B"  represents  quality 
slightly  below  standard,  etc.  to  "P  representing  material  not  likely  to  survive  a  poling 
operation.  "Scatter"  is  an  estimate  of  the  number  of  particles  in  a  25  mm  length  of 
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7a/)fe  3.  Summary  of  growth  runs. 


crystal  detectable  with  the  unaided  eye  using  a  white  microscope  light  in  a  darkened 
room. 

Macroscopic  crystal  quality  was  found  to  be  a  strong  ftinction  of  the  growth  axis.  Z*axis 
material  was  the  easiest  to  grow  with  low  strain.  Hie  best  boules  (B  and  B+  on  our 
strict  grading  scale)  were  only  slightly  inferior  to  undoped  LiNbOj  used  in  Q-switch 
applications.  With  y-axis  and  especially  x-axis  growth,  it  was  discovered  that  higher 
effective  thermal  gradients  were  necessary  to  avoid  the  formation  of  polycrystaliinity.  As 
a  consequence,  the  resulting  boules  tended  to  be  hetivily  strained.  All  three  x-axis  boules 
were  cracked,  two  of  them  severely  enough  to  preclude  poling. 

In  general,  for  a  given  axis  of  growth,  a  higher  MgO  melt  concentration  resulted  in  lower 
macroscopic  crystal  quality.  No  effect  was  observed  of  the  relatively  small  variation  in 
Li/Nb  melt  ratio. 

All  of  the  grown  crystals  exhibited  planar  striae  perpendicular  to  the  growth  axis  similar 
to,  although  less  severe  than,  those  described  in  the  Nightingale  work  It  is  notable 
that  the  striae  were  only  visible  in  light  with  a  polarization  component  along  the 
crystallographic  z-axis.  For  light  polarized  in  the  crystallographic  x-y  plane,  the  striae 
were  invisible.  Heat  treatment  up  to  50  hours  in  duration  at  temperatures  in  excess  of 
1200®C  appeared  to  have  no  influence  on  reducing  the  striae.  The  appearance  of  the 
striae  did  not  vary  strongly  with  MgO  concentration  over  the  range  studied. 

Chemical  analysis  (subcontracted  to  Ledoux  &  Co.  of  Teaneck,  NJ)  of  several  boules 
indicated  an  effective  distribution  coefficient  for  Mg  of  roughly  1.05.  Thus,  MgO 
concentration  decreased  down  the  growth  axis  of  each  boule. 

IR  transmission  analysis  produced  surprising  results.  As  expected  from  Bryan’s  report’, 
boules  with  the  lowest  MgO  concentration  exhibited  unshifted  (2.87m)  OH'  absorption 
peaks  while  boules  with  the  highest  MgO  concentration  exhibited  shifted  (2.83m)  peaks. 
For  boules  with  intermediate  MgO  concentration,  however,  boule  bottoms  tended  to 
have  shifted  OH'  absorption  peaks  while  boule  tops  tended  to  have  unshifted  peaks. 
Given  the  results  of  the  chemical  analysis  indicating  a  decrease  in  MgO  concentration 
down  the  length  of  a  boule,  these  findings  suggest  that  the  OH'  absorption  peak  shifting 
phenomenon  is  not  a  simple  function  solely  of  MgO  concentration. 


3.1.1  Ferroelectric  poling 

Ferroelectric  poling  of  LiNbOj  is  normally  accomplished  by  imposing  an  electric  field 
along  a  crystal’s  polar  z-axis  while  cooling  it  thiough  its  Curie  transition.  The  poling 
operation  can  be  incorporated  into  the  Czochralski  proce.ss  itself  or  dealt  with  in  a 
separate  post-growth  proce.ssing  step.  The  effectiveness  of  a  given  set  of  poling 
conditions  is  most  easily  tested  through  standard  techniques  for  LiNb03  domain 


decoration  involving  polishing,  etching,  and  optical  microscopy®.  Undoped  LiNbOj  of  the 
congruent  composition  has  a  Curie  temperature  of  roughly  1140°C  and  a  melting  point  of 
roughly  1250“C.  Typical  poling  conditions  involve  a  poling  field  of  roughly  0.5  V/cm. 

In  the  present  work,  poling  was  carried  out  in  a  separate,  post-growth  process.  Actual 
poling  experimentation  was  limited  to  the  first  two  crystals  grown.  As  shown  in  Table  3, 
Curie  temperature  values  for  the  MgOrliNbOj  grown  in  this  work  ranged  from  1205  to 
1222‘*C.  Given  an  MgOrUNbOj  melting  point  less  than  that  for  undoped  LiNbOj,  the 
range  of  poling  process  starting  temperature  is  comparatively  limited. 

Despite  the  consequent  constraint  of  necessarily  well-controlled  temperature 
programming  in  the  poling  furnace,  it  was  determined  that,  provided  the  entire  boule 
was  allov'ed  sufficient  time  to  reach  a  temperature  throughout  its  volume  in  excess  of  its 
Curie  temperature,  the  0,5  V/cm  poling  field  used  for  undoped  LiNbOj  was  sufficient  to 
pole  \he  MgOiUNbOj  without  attendant  damage.  Although  a  number  of  boules  were 
checked  using  polishing/etching/microscopy  for  the  effectiveness  of  the  poling  operation, 
there  was  no  evidence  of  domain  multiplicity  either  in  these  tests  or  in  subsequent 
optical  characterization  of  any  of  the  crystals. 


3.1,2  Photorefractive  effect  (PRE)  testing 

A  number  of  the  crystals  were  characterized  for  PRE  susceptibility.  Samples  from  Z-6, 
Z-8,  Y-2,  Y-3,  and  a  standard  undoped  optical-grade  LiNbOj  boule  were  fabricated  or 
re-fabricated  into  test  specimens  measuring  5.5x7x8mm  (x  x  y  x  z).  The  x-faces  were 
polished  to  typical  optical  surface  specificatiora  but  not  anti-reflection  coated.  The  test 
set-up  is  shown  in  figure  22.  An  argon  ion  laser  with  3.2W  of  cw  operating  power  at 
488nra  was  used  in  creating  PRE  damage  in  the  specimens.  Beam  propagation  was 
along  X  with  polarization  parallel  to  z.  PRE  damage  was  quantified  through 
measurement  of  optical  power  passing  through  an  aperture  placed  after  the  crystal  in  the 
optical  path.  For  each  specimen,  PRE  damage  was  measured  versus  time,  temperature, 
and  incident  power. 

Figure  23  illustrates  the  results.  Each  data  point  represents  the  normalized  laser  power 
passing  through  the  aperture  after  reaching  steady-state  PRE  conditions  in  a  given 
sample.  After  a  change  in  temperature  or  laser  fluence,  tens  of  seconds  were  typically 
required  to  reach  such  steady-state  conditions.  No  hysteresis  was  detected  on  cycling  the 
temperature  or  laser  fluence,  provided  sufficient  time  was  allowed  to  reach  steady-state. 

Several  points  are  notable  about  the  data  in  figure  23.  First,  for  a  given  sample  and  a 
given  fluence,  thermal  annealing  of  PRE  is  evident.  Contrary  to  previous  reports, 
however,  the  "annealing  temperature"  ensuring  absence  of  PRE  is  clearly  a  function  of 
both  MgO  concentration  and  fluence.  Second,  PRE  in  a  given  sample  at  a  given 
temperature  is  a  function  of  laser  fluence.  Although  not  surprising,  this  observation  has 
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Figure  21  Photorefractive  effect  (PRE)  testing. 
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Tho  Maximum  Inserted  Fluence  Is  43  Kliowaffs  Per  Square  Centimeter.  PRE  Causes  Increased 

Beam  Divergence  (Etpecialiy  Parallel  To  The  Z-axIs)  The  Result  1$  Increased  Loss  At  Aperture 
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Figure  23.  PRE  in  Mg.LiNbOj  vs.  temperature,  fMgO]  and  beam  intensity. 
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F^ure  24.  Phasematching  temperature  as  a  function  of  fMgOJ. 
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not  been  previously  reported.  Tliird,  similar  to  Bryan’s  report  at  a  given  temperature 
and  fluence,  PRE  is  clearly  a  function  of  MgO  concentration  for  concentration  values 
less  than  roughly  5  mol%.  Although  the  data  were  left  off  figure  23,  the  test  specimens 
fabricated  from  undoped  UNbOj  were  observed  to  allov/  minimum  power  through  the 
aperture  for  a  43  kW/cm^  fluence  until  the  first  onset  of  thermal  annealing  effects  above 
IbS^C.  Fourth,  for  the  one  sample  tested  with  a  sufficiently  high  MgO  concentration  to 
exhibit  a  "shifted"  IR  absorption  peak,  it  was  impossible  to  induce  measurable  PRE  for 
all  tested  conditions  of  fluence,  exposure  time,  and  temperature.  The  presumably  most 
severe  set  of  these  conditions  involved  a  fluence  of  43  kW/cm^  for  16  hours  at  20*C. 


3.13  NLO  phase-matching  measurements 

Measuremente  were  made  of  the  temperature  for  non-critical  phase-matching  of  second 
harmonic  generation  from  a  1064  nm  laser  source.  The  subject  crystals  were  the  same  as 
those  characterized  for  PRE  susceptibility  and  the  test  specimens  were  prepared  to  the 
same  dimensions  and  specifications.  A  single  longitudinal  mode  laser  was  Q-switched  at 
10  Hz  to  produce  temporally  smooth  pulses  of  11  ns  duration.  Propagation  was  along  x 
with  the  1064  mn  fundamental  polarized  along  z  and  the  532  nm  second  harmonic 
polarized  along  y.  The  fundamental  and  second  harmonic  were  separated  using  a 
thin-film  interference  filter.  Signal  intensity  at  each  wavelength  was  measured  with  fast 
silicon  detectors.  Sample  temperature  was  controlled  and  ramped  in  cycles  using  a  small 
oven  and  controller.  Sample  temperature  was  measured  with  a  type-J  thermocouple. 
Phase-matching  temperature  was  determined  as  the  ratio  of  532  nm  signal  to  1064  nm 
reached  its  highest  value  in  each  temperature  cycle. 

The  results  are  shown  in  figure  24.  The  dramatic  peak  in  phase-matching  temperature 
versus  MgO  concentration  has  not  previously  been  reported.  It  is  notable  that  the  MgO 
concentration  values  in  the  region  of  the  peak  correspond  to  the  discontinuous  shift  in 
the  IR  absorption  peak  for  OH‘.  A  physical  explanation  for  such  behavior  is  not 
obvious,  although  similar  peaks  in  the  composition  dependence  of  phase-matching 
temperature  have  recently  been  reported  for  vapor  transport  equilibrated  LiNbOj 
crystals  with  solid  compositions  near  50  mol%  Li20’.  These  of^congruent  composition 
crystals  were  not  tested  for  photorefractive  effects. 
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In  terms  of  macroscopic  crystal  quality,  MgOrLiNbOj  is  clearly  more  difficult  to  grow 
than  undoped  LiNbOj,  and  the  difficulty  increases  with  MgO  concentration.  For  a  given 
MgO  concentration,  growth  along  the  crystallographic  z-axis  is  easier  than  growth  along 
the  y-axis  which  is  in  turn  easier  than  growth  along  the  x-axis.  Planar  striae  are  present 
for  all  growth  orientations.  The  striae,  with  a  spacing  of  roughly  200  pirn,  run 
perpendicular  to  the  growth  direction  and  are  visible  with  z-polarized  light. 


The  ferroelectric  poling  conditions  that  were  developed  for  Mg0;LiNb03  based  on 
standard  poling  conditions  for  undoped  LiNbOs  adjusted  for  the  increased  Curie 
temperature  are  effective  and  reproducible. 

PRE  damage  threshold  is  a  function  of  temperature  and  laser  fluence.  Both  PRE 
damage  threshold  and  NLO  phase-matching  temperature  are  fiinctions  of  MgO 
concentration,  particularly  in  the  vicinity  of  5  mo\%.  For  one  crystal  doped  with  a  high 
concentration,  it  was  impossible  to  induce  measurable  PRE  damage  for  all  tested 
conditions  of  laser  power,  exposure  time,  and  temperature. 

The  discontinuous  shift  in  the  IR  absorption  peak  associated  with  residual  OH'  appears 
to  be  an  important  indicator  of  some  dramatic  optical  property  changes  in  Mg0:LiNb03. 
The  peak  shift  appears  to  be  a  function  of  several  parameters  including  but  not  limited 
to  MgO  concentration. 

Commercial  NLO  product  development  efforts  are  currently  underway  utilizing  the  PRE 
damage-resistant  MgOiLiNbOj  composition  identified  through  this  projea.  Of  particular 
interest  is  resonant  frequency  doubling  of  a  1064  nm  fundamental.  Various 
waveguide-based  NLO  technologies  are  under  development  as  well. 

Additional  research  on  Mg0:LiNb03  is  warranted  to  investigate  several  issues  identified 
in  this  project.  First,  the  relationships  between  crystal  composition,  point  defect 
structures,  and  PRE  mechanisms  appear  to  be  both  theoretically  complex  and 
technologically  significant.  Second,  the  effect  of  MgO  additions  beyond  the  5  mol% 
level  remains  largely  unexplored.  Third,  if  desirable  from  the  standpoint  of  optical 
performance  (even  though  scattering  levels  are  low  enough  for  present  applications),  the 
possibility  exists  for  further  reduction  or  even  elimination  of  the  observed  planar  striae 
through  application  of  Czochralski  process-control  improvements  currently  under 
development  for  undoped  LiNb03. 

An  issue  not  explored  in  this  work  concerns  the  influence  of  electric  fields  on  PRE 
damage  in  MgOrLiNbOj.  A  significant  decrease  in  PRE  damage  threshold  has  been 
observed  in  some  Mg0:LiNb03  subjected  to  DC  fields  for  purposes  of  piezoelectric 
locking  to  laser  cavity  resonances.  A  complete  understanding  of  PRE  in  MgOrUNbOs 
will  have  to  account  for  such  observations. 

Finally,  it  should  be  mentioned  that  work  is  in  progress  doping  LiNb03  with  other 
species  analogous  to  MgO'°.  Through  this  work,  it  is  hoped  to  maintain  the  beneficial 
effects  of  PRE  resistance  while  avoiding  or  reducing  some  of  the  processing-related 
deficiencies  of  Mg0:LiNb03. 
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CONCLUSIONS 


It  appears  from  the  present  work  that  absorption  losses  are  significantly  higher  than 
scattering  losses  in  Mg0:UNb03,  regardless  of  MgO  concentration,  as  long  as  no 
catastrophic  flaws  (cracks)  are  visible.  The  likelihood  of  flaws  increases  with  MgO 
concentration,  although  X-grown  crystals  are  the  most  likely  to  present  flaws,  even  with 
small  amounts  of  MgO.  The  quality  of  Z-grown  boules  appeared  more  immune  to 
[MgOj,  and  Y-grown  boules  showed  sensitivity  to  [MgO]  between  that  of  X-  and  Z* 
grown  material. 

The  changes  in  formulation  in  the  intrinsic  crystal,  where  (Li2^]  is  changed  between 
48.60%  to  4838%  appear  to  have  no  effect  in  absorption,  or  at  least  are  drowned  by 
variations  observed  in  crystal  quality.  It  appears  then  that  this  slight  change  in 
formulation  will  not  affect  other  conclusions  that  can  be  drawn  from  the  data. 

The  optimal  results  in  scattering  and  absorption  were  obtained  with  Z-grown  crystals, 
perhaps  due  to  prior  experience. 

The  most  striking  observation  concerns  the  means  of  attaining  low  absorption:  using  Z- 
polarized  light,  whereas  Z-propagation  causes  the  most  absorption.  This  should  have 
important  consequences  for  the  applications  as  well  as  the  choice  of  growth  axis, 
depending  on  the  total  light  path  required. 

Finally,  the  fact  that  large  variations  can  be  observed  among  different  samples  from  the 
same  boule  points  to  the  lack  of  unifomiity  found  in  this  kind  of  crystal.  In  most  cases, 
the  variations  within  one  boule  are  comparable  to  the  differences  among  boule  averages 
-both  in  scattering  and  absorption.  In  cases  where  the  quality  is  critical,  measuring 
attenuation  in  each  crystal  may  be  the  only  way  to  insure  reliability. 
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APPENDIX  A  -  THERMAL  ANALYSIS 

Consider  light  incident  on  a  crystal,  with  incident  power  Pq.  The  injected  power  is  then 
Pj  =  T  (Al) 


where  T  »  1-R  is  the  transmission  coeFicient  across  the  surface,  and  the  reflecti\aty  is 
givrn  by  the  index  of  refraction  of  the  crystal,  assuming  normal  incidence: 

R  =  ((n  - 1)  /  (n  T  1)1“  (A2) 

after  traveling  through  a  length  of  crystal  L,  the  absorbed  power  is 

A,  =TP,(l-e-“)  (A3) 

and  upon  reaching  the  other  end  of  the  crystal  a  fracti'.  .  ‘.he  remaining  power  will  be 
reflected  back,  yielding  a  second  contribution  to  absorption: 

A,  =  T  P,  e  “  R  (  1 .  e  '  ■• )  (A4) 
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Adding  all  the  contributions  from  internal  reflections  leads  to  a  total  absorbed  power 


A  =  Aq  +  Aj  +  ...  »  K  Pq  (A5) 

where  k  is  the  absorbed  fraction,  and  the  absorption  coefficient  X  can  be  expressed  in 
terms  of  k,  T  and  R  by 

X  =(l/L)ln{(T.KR)/(T.K)}  (A6) 

which  in  the  T  «  1,  R  «  0  approximation  reduces  to 

A  = -(1/L)ln(l -k)  (A7) 


Calorimetry 

Now  let  us  consider  a  body  of  mass  m  and  heat  capacity  Ce-  The  temperature  change 
dT/dt  is  given  by  its  change  in  heat: 

dq/dt  =  m  Ce  (dT/dt)  (A8) 

which  is  equal  to  the  fraction  k  of  absorbed  power  Pq  minus  the  radiated  power: 

m  Ce  dT/dt  =  k  Po  -  p  (T-To)  (A9) 

where  p  is  a  radiation  coefficient  accounting  for  radiative  and  conductive  thermal  losses, 
and  Tq  is  the  ambient  temperature.  Changing  variables  to  relative  temperature,  the 
equivalent  of  subtracting  the  ambient  temperature  as  monitored  by  a  neighboring 
thermistor, 


y  =  t-To 

(AlO) 

equation  (9)  above  can  be  written  as 

dy/dt  =  a  Po  -  r  y 

(All) 

with  an  absorption  term  a  =  k  /  (m  Cg)  and  a  radiative  term  r  =  p  /  (m  Cg).  If  one 
heats  the  body,  the  absorption  and  radiative  terms  are  present,  yielding  a  time 
dependency  (defining  y(0)=0): 


y(t)  =  Po  (a/r)  (  1  -  e"') 


(A12) 


but  during  cooling  only  the  radiative  term  is  present,  so  if  at  t’=0  we  define  y(t’)®yot 


y(t’)  =yoe'^'’  (A13) 

In  other  words,  the  term  r  can  be  obtained  from  the  cooling  off  curve,  and  then  used  to 
determine  the  term  a  from  the  heating  curve.  The  cooling  off  curve  is  first  linearized 
and  then  fit  to  a  straight  line  by  a  least  squares  algorithm:  the  slope  determines  r  since 

In  (  y(t’)/yQ )  =  -  r  f  (A14) 

then  the  heating  curve  can  be  Uuearized  and  least-squares  fit  to  obtain  a: 

y(t)r/Po  =  a(l-e-^‘)  (A15) 

This  allows  taking  into  account  all  data  points,  and  a  judicious  sampling  allows  a 
reasonable  weighing  of  the  samples. 

Thermistor  mass  corrections 

To  account  for  the  contribution  of  the  thermistor /epoxy  mass  m,  to  the  thermal  mass  of 
the  crystal  m,.,  we  can  define  the  equivalent  mass 

m,q  =  m,  +  m,  X  (Ce,  /  C^) 

where  Cg,  is  the  thermistor/epoxy  heat  capacity  and  Ce^,  that  of  the  crystal  by  itself.  All 
the  crystals  were  weighed  before  and  after  applying  the  thermistor,  so  ra,  is  known.  Ce, 
was  estimated  at  10.6  (Cal/®-mol),  the  value  for  Si02,  as  recommended  by  the  thermistor 
manufacturer. 
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APPENDIX  B  -  FRESNEL  EFFECTS 


When  a  beam  of  light  enters  normal  to  the  surface  of  a  crystal  whose  exit  surface  is 
nearly  parallel  to  the  first  one,  a  set  of  Fresnel  reflections  is  established.  Since  for  the 
case  at  hand  the  reflection  coefficients  are  rather  small,  the  interference  pattern  is  of  no 
concern. 

Consider  an  incident  beam  with  intensity  striking  a  crystal  of  thickness  d,  whose 
interface  has  a  reflectivity  R  (transmissivity  T  =  1  -  R  ),  as  shown  in  figure  B-1.  The 


Figure  B~l.  Reflected  beam  series  in  a  cell  with  parallel  walls. 
reflected  beam  intensity  can  be  expressed  as  a  series  that  adds  up  to 


1,  =  I„  R  {  1  +  (  1  -  2  R)  }  /  (1  -  RV]  (B1) 

where  L  =  e'^  (B2) 

and  X  is  the  total  loss  (absorption  a  and  scattering  a): 

X  =  a  +  o  (B3) 

Similarly,  the  transmitted  intensity  is  a  series  that  converges  to 
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It  =  lo  L  /  {1  -  RV] 


(B4) 


and  the  lost  (absorbed  or  scattered)  intensity  (after  infinite  number  of  internal 
reflections)  is 


Ia  =  IoT(1-L)/(1-LR)  (B5) 

which  is  the  equivalent  of  equation  A-6  in  Appendix  A. 

Let  us  now  look  a  the  beam  intensity  as  it  zig-zags  through  the  crystal.  The  remaining 
intensity  after  travelling  a  distance  /  =  n  d  +  x  is 

1(0  =  lo  T  R"  e^  (B6) 

that  is,  it  exponentially  decays  along  the  material,  suffering  drastic  losses  at  each 
reflection.  This  dependency  is  shown  in  figure  B-2,  for  a  rather  lossy  material  (a  =  0.5) 
one  inch  thick. 
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B-2  Remaining  beam  intensity. 


The  fraction  of  the  beam  intensity  that  is  absorbed  at  the  interval  between  x  and  x+dx, 
(provided  the  interval  does  not  contain  a  reflexion),  is 

l3(/, /+dx)  =  1(0  X  a  X  dx  (B7) 
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Absorbed  frocfion 

l.«nglh==25  mm,  Lombdo^O.AST/m,  n=2.25 


'  R^fTusInlno  intensity  Absorbed  power 


Figure  B-3.  Absorbed  fractiorL 

Since  1(0  is  discontinuous,  the  integral  is  a  rather  un-elegant  series.  It  was  calculated  for 
the  same  case  described  above,  and  it  is  shown  in  figure  B-3.  Note  that  figures  B-2  and 
B-3  contain  the  same  curves,  differing  only  in  the  venical  scale. 

Although  the  intensity  of  the  beam  decays  rapidly  within  a  few  reflections,  it  remains 
large  compared  with  the  absorbed  fraction:  if  the  beam  touches  the  walls  generating 
absorption  against  non-transparent  walls,  even  after  four  reflertions,  this  contribution  to 
the  deposited  power  can  cause  significant  measurement  errors.  This  effect  was  an 
important  consideration  in  performing  the  measurements  outlined  in  the  present  report. 
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